Time-odd mean fields in the rotating frame: 
microscopic nature of nuclear magnetism. 

A. V. Afanasjevi and P. Ring 
Physik- Department der Technischen Universitdt Miinchen, D-85747 Garching, Germany 

(February 8, 2008) 



Abstract 



The microscopic role of nuclear magnetism in rotating frame is investigated 
for the first time in the framework of the cranked relativistic mean field theory. 
It is shown that nuclear magnetism modifies the expectation values of single- 
particle spin, orbital and total angular momenta along the rotational axis 
effectively creating additional angular momentum. This effect leads to the 
increase of kinematic and dynamic moments of inertia at given rotational 
frequency and has an impact on effective alignments. 
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It is a well established fact that microscopic mean field theories provide a good description 
of the properties of finite nuclei. The most succesful among them are the non-relativistic 
theories based on the finite range forces of the Gogny type and zero-range forces of the 
Skyrme type and relativistic mean field (RMF) theory 0. In these theories the physical 
observables are most sensitive to the time-even fields, which has widely been investigated 
throughout the periodic table. On the other hand, the properties of time-odd fields, which 
appear only in the nuclear systems with broken time-reversal symmetry, are poorly known. 
However, it is known nowadays that these fields are very important for the proper description 
of rotating nuclei |3|-§|], magnetic moments J7|, the properties of the N = Z nuclei |§ and 
pairing correlations ||. 

In rotating nuclei, the time-odd mean fields emerge both from the Coriolis operator and 
from the parts of the Hamiltonian (Lagrangian) related to the currents. Both in relativistic 
and non-relativistic ||,|5,|10| mean field theories the most important modification coming 



from the presence of currents lies in the renormalization of the moments of inertia. They 
increase the kinematic (J^) and dynamic (J^) moments of inertia by ~ 20—30%. However, 
the microscopic mechanism of these modifications has not been well understood till now. 

The aim of the present article is to investigate this mechanism staying strictly in the 
framework of the cranking model. As a theoretical tool we are using the cranked relativistic 
mean field (CRMF) theory HIB in wnich the time-odd mean fields are defined in a unique 



way. Since we are interested in a general understanding of the microscopic mechanism of 
the renormalization of the moments of inertia, the pairing correlations are neglected in the 
present investigation. 

In RMF theory the nucleus is described as a system of point-like nucleons (Dirac spinors) 
coupled to the mesons and to the photons. The nucleons interact by the exchange of several 
mesons: isoscalar scalar cr-meson, isoscalar vector w-meson, isovector vector p-meson and 
the photon. The transformation to the rotating frame in the one-dimensional cranking 
approximation leads to the CRMF theory, the details of which are given in Refs. |TT| , |I2"| ,P||S|] . 
The main features of this theory important in the present discussion are outlined below. 



In the Hartree approximation, the stationary Dirac equation for the nucleons in the 
rotating frame is given by 

| a (_, V - V(r)) + V (r) + (5{m + S(r)) - Sl x j x } A = 

= €i1pi (1) 

where Vo(r) represents a repulsive vector potential: 

1 -T 3 

V (r) = g w u {r) + g p T 3 p (r) + e - A (r) (2) 

which contains time-like components of the vector mesons and S(r) an attractive scalar 
potential: 

S(r) = g a a(r). (3) 

A magnetic potential V(r): 

V(r) = g uU ( r ) + g p T 3 p(r) + e^p-A(r) (4) 

originates from the space-like components of the vector mesons. Note that in these equations, 
the four- vector components of the vector fields p M and N 1 are separated into the time-like 
(uq, po and A ) and space-like [u> = (u; x , uj z ), p = (p x ,p y ,p z ), and A = (A 21 , A 2 )] 
components. Finally the term 

~~ Qx-Jx — —Qx(L x + t^ 11 ) (^) 

represents the Coriolis field. 

The time-independent inhomogeneous Klein-Gordon equations for the mesonic fields are 
given by 



{-A - (n x L x f + ml} a(r) = -g c [tf (r) + p n s {r) 



-g 2 cr 2 (r) - g 3 a 3 (r) 



{-A - (Q X L X ) 2 + ml} wb(r) = [tf(r) + ^(r)] , 
{-A - (Q X (L X + S x )) 2 + ml} «(r) = </ w tf*(r) + j»] , 



(6) 



with source terms involving the various nucleonic densities and currents 

pr(r) = f>*W) f AMr), 
i=i 

= E(V'i(r)) t V« i (r) I 
i=i 

j»*(r) = X;(^(r)) t a^ i (r), (7) 
i=i 

where the labels n and p are used for neutrons and protons, respectively. In the equations 
above, the sums run over the occupied positive-energy shell model states only (no-sea ap- 
proximation) [13|]. For simplicity, the equations for the p meson and the Coulomb fields are 
omitted in Eqs. (|6]) since they have the structure similar to the equations for the u meson. 
Since the coupling constant of the electromagnetic interaction is small compared with the 
coupling constants of the meson fields, the Coriolis term for the Coulomb potential Ao(r) 
and the spatial components of the vector potential A(r) are neglected in the calculations. 
Note that in the CRMF theory the currents j n,p (r) are the products of the large (G) and 
small (F) components of the Dirac spinors ipi(r). 

The Coriolis operator J x and the magnetic potential V(r) in the Dirac equation as well 
as the currents j n ' p {r) in the Klein-Gordon equations do not appear in the RMF equations 



for time-reversal systems |Tj] since they break this symmetry [fLq| . Similar to the non- 
relativistic case Q, the presence of the Coriolis operator leads to the appearance of the 
time-odd mean fields. However, the CRMF calculations with only these time-odd fields 
accounted underestimate the experimental moments of inertia by 20-30% @|J. A similar 
situation holds also in non-relativistic theories, see Refs. [§,|1(J an d references therein. The 



inclusion of the currents j n ' p (r), which leads to the space-like components of the vector u 
and p mesons and thus to magnetic potential V(r), considerably improves the description 
of experimental moments of inertia. The effect coming from the space-like components of 
the vector mesons is commonly referred as nuclear magnetism (further NM) |Tl| in the 
framework of RMF theory. 

Before we discuss the microscopic impact of NM on the moments of inertia, we remind 



that the kinematic (J^- 1 ') and dynamic (J*- 2 -*) moments of inertia are defined as 

J(1)( ^ ) = £' J(2)( ^ ) = ^' (8) 

where the rotational frequency Q x along the x-axis in a one-dimensional cranking approx- 
imation is determined by the condition that the expectation value of the total angular 



momentum J at spin I has a definite value [16 



J=($n\ J, I $n> = v // ( / + !) = £(* \j* V ( 9 ) 

i=i 

Note that J is defined in the cranking model as a sum of the expectation values of the 
single-particle angular momentum operators j x of the occupied states. This suggests that 
the modifications of the moments of inertia due to NM should be traced to the changes of 
the single-particle (j x )i values. 

In order to confirm this expectation the CRMF calculations with and without NM (the 
later will be further denoted as WNM) have been performed for the doubly magic superde- 
formed (SD) configuration in 152 Dy (7r6 4 z/7 2 ) ||. This configuration has been selected due 
to the large stability of the SD minimum against rotation. The calculated (j x )i values of 
the single-neutron orbitals forming the N = 1 shell are given in Fig. |I|. This shell has been 
selected because (i) it is reasonably well separated in energy from the N = 3 shell and thus 
the emerging picture is not strongly disturbed by an interaction between the N = 1 and 
the N = 3 shells, (ii) the number of single-particle orbitals is, from one hand, reasonably 
small to make an analysis transparent and, from other hand, sufficient to make meaningful 
conclusions. 

Indeed NM has a considerable impact on the expectation values of the single-particle 
angular momentum (j x )i- If we define the contribution to (j x )i due to NM as 

AO,), = &>f M - Q X )T NM (io) 

then one can see in Fig. [I] that A(j x )i is positive at the bottom and negative at the top of 
the N = 1 shell. The absolute value of A(j x ) correlates with the absolute value of (j x )i- 



Note that the contributions to (j x )i due to NM are small in the middle of the shell. It was 
checked also that similar features hold also in higher- N shells. Note that the total angular 
momentum built within the N = 1 shell is not equal (see top panel in Fig. P due to the 
admixture of the N = 3 shell and that NM increases mixing of the N = 3 and N — 1 shells 
as seen from an increased value of this momentum and from the analysis of the structure of 
the wave functions of the single-particle orbitals. These results suggest that if it would be 
possible to isolate the iV-shell then the contribution of NM to total angular momentum of 
this shell would be equal zero. 

In order to gain a deeper understanding of the impact of NM, the angular momentum 
of the particle along the rotation axis is separated into the orbital and the spin parts: 
jx = 4 + s x . Then the contributions due to NM to the expectation values of the orbital 
(A{l x )i) and the spin (A(s x )i) angular momenta, which are defined similar to Eq. (|T0|), are 
plotted in Fig. ||. 

The spin angular momenta of the particles in the [110]l/2 + and [110] orbitals, which 
are the lowest orbitals in the N — 1 shell, are almost fully aligned along the axis of rotation 
already at fl x = 0.0 MeV. Their expectation values (s x )i are equal to ~ 0.49ft and ~ — 0.49ft 
at all rotational frequencies employed. Thus NM cannot have an impact on the spin angular 
momenta of the particles in these orbitals (A(s x )i ~ Oft, see Figs. ||a,b) and the contributions 
to (j x )i due to NM are almost solely defined by the modifications of the expectation values 
of the orbital angular momentum (l x )i induced by NM, see Figs. ^|a,b. On the contrary, 
the spin angular momenta of the particles in the other orbitals forming the N — 1 shell 
are only partially aligned along the rotation axis. Thus NM modifies the expectation values 
of both the spin and orbital angular momenta of the particles in these orbitals, see Figs. 
|2]c,d,e,f. Note that the contributions A(l x )i and A(s x )i have in general different rotational 
frequency dependence and can have different signs. The later is clearly seen in the case of 
the [101]3/2 + orbital, see Figs. ^|c. The results for proton orbitals are very similar to the 
ones presented here for neutrons. Thus one can conclude that, in general, NM leads to the 
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modification of alignment properties of single-particle orbitals on the level of spin, orbital 
and total angular momenta. 

It was checked that the impact of NM on (j x )i has the same general features also in other 
shells. However, for higher iV-shells the crossings between the single-particle routhians make 
the emerging picture more complicate. The question is then how to understand the increase 
of kinematic and dynamic moments of inertia induced by NM. In the case of yrast SD band 
in 152 Dy, the shells up to N = 3 can be considered as 'closed' since all the orbitals belonging 
to them are occupied. The contribution of the 'closed' shells to the expectation value J of the 
total angular momentum accounts only for few % of J as illustrated on the example of the 
N = 1 shell, see top panel in Fig. [j]. Thus the main contribution to J (and correspondingly 
to j( l > and J^ 2 ') is coming from the particles and holes outside the 'closed' shells. 

In the lowest SD configuration of 152 Dy, valence protons are sitting in the orbitals, emerg- 
ing from the bottom of the N = 5 and N = 6 shells (see for example Fig. 2 in Ref. Jl7|]), for 
which the contribution of NM to (j x )i is positive. Valence proton holes are located in the 
orbitals, emerging from the middle and the top of the N = 4 shell. The contribution of NM 
to (j x )i for these orbitals is either close to zero or negative. However, these contributions 
should be taken with opposite sign due to hole nature. Thus these proton holes also give 
a positive contribution due to NM to the total proton angular momentum. The consider- 
ation of the neutron subsystem can be performed in a similar way and the conclusions are 
the same. Thus one can conclude that at given rotational frequency the modifications of 
single-particle (j x )i values induced by NM will lead to an increase of total J and, as a result, 
to an increase of kinematic and dynamic moments of inertia. Essentially this means that 
NM is the source of the creation of angular momentum in rotating nuclei in addition to the 
Coriolis force. Clearly this increase in J, and induced by NM should depend both 
on the rotational frequency and on the shell filling (and thus on the specific single-particle 
configuration). 

In the light of these findings, it is important to see which other physical quantities are 
affected by NM. NM can have considerable impact on the signature splitting of the signature 
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partner orbitals as it is clearly seen in the case of the [110]l/2 ± and [330] l/2 ± orbitals, see 
Fig. [| One should also note that the change in (j x )i induced by NM is not fully reflected in 
the change of the slope of the single-particle energies drawn as a function of the rotational 
frequency (see Fig. This is contrary to the results of the phenomenological models based 
on the Woods-Saxon or Nilsson potentials where for a cranking at a fixed deformation, 
the slope of the orbital routhian against rotational frequency corresponds to the alignment 
(jx)i — —Se i x /8Q X , where e i x is the single-particle energy of i-th orbital at given Q x . For 
example, the [110]l/2 ± and [101]l/2~ orbitals have the contributions to induced by NM 
(see Figs. ||a,b,f) similar in absolute value. However, the change of the slope of the single- 
particle routhians induced by NM is larger for the [110]l/2 + orbital as compared with the 
[110]l/2 + and [101]l/2~ orbitals, see left panel of Fig. |[ This essentially suggests that the 
difference between the single-particle energies of the orbitals calculated with and without 
NM reflects not only the change of alignment properties of corresponding single-particle 
orbitals but also the energy gain (or cost) due to the interaction of the nucleon with the 
magnetic potential V{r). 

In order to check that the features mentioned above are not coming from small differences 
in equilibrium deformations obtained in the calculations with and without NM, the calcu- 
lations without NM have been performed starting from the self-consistent fields obtained 
in the calculations with NM and performing only one iteration during which the physical 
quantities of interest have been recalculated neglecting NM. Although the solution obtained 
in this case is not self-consistent, its equilibrium deformation differs only marginally from the 
one obtained in fully self-consistent calculations with NM. Using this result a similar analysis 
as presented above has been repeated. Although some minor quantitative differences exist, 
the physical conclusions are the same. 

It is reasonable to expect that also the physical observables which depend on the align- 
ment properties of specific single-particle orbitals are affected by NM. One of these ob- 
servables is the effective alignment i e ff which measures the difference between the spins 
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of bands A and B at constant Q x : i e ^(fi, x ) = Ib(Q x ) ~~ Ia(&x) and which in Cranked- 
Nilsson+Strutinsky (CNS) approach based on the cranked Nilsson potential is predomi- 
nantly defined by the alignment (s) of the single-particle orbital (s) by which the compared 
bands differ |H| . Fig. |] shows that NM indeed can have a considerable impact (especially 
in the case of the 7r[651]3/2(r = +i) orbital) on effective alignments. A similar conclusion 
about the impact of time-odd components of the mean field on relative (effective) alignments 
has been drawn earlier in the cranked Hartree-Fock approach based on Skyrme forces ||. 
Note that contrary to the CNS approach polarization effects play a much more important 
role in the CRMF theory and thus effective alignments are not defined predominantly by 
single-particle alignments, see Fig. 01. The fact that effective alignments between SD bands 



in different mass regions []1Q| plfl are correctly described in the CRMF theory indicates that 
both alignment properties of single-particle orbitals and relevant polarization effects are 
properly accounted in this theory. 

In conclusion, the microscopic role of nuclear magnetism in rotating frame has been in- 
vestigated for the first time. The breaking of time-reversal symmetry induced by rotation 
results in baryonic currents which lead to the non-vanishing space-like components of the 
vector mesons creating the magnetic potential entering into the Dirac equation. The mag- 
netic field created by this potential modifies the expectation values of the single-particle 
angular momenta effectively creating additional total angular momentum. Both the spin 
and the orbital single-particle angular momenta are affected by nuclear magnetism. This 
effect is state- and rotational frequency dependent. It leads to the increase of total angular 
momentum, kinematic and dynamic moments of inertia at given rotational frequency and 
has an impact on signature splitting, single-particle energies and on effective alignments. 
It is reasonable to expect that an increase in kinematic and dynamic moments of inertia 
resulting from the time-odd components of the mean field in the non-relativistic mean field 
theories flJTOfl is caused by a similar microscopic mechanism as the one discussed in the 
present manuscript. Clearly more detailed studies of the microscopic nature of nuclear mag- 
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netism both in non-rotating and rotating cases are mandatory. Such an investigation is in 
progress and its results will be reported later. 
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FIG. 1. The expectation values (j x )i of the neutron orbitals forming the N = 1 shell calculated 
with and without NM. They are given along the deformation path of the lowest SD configuration 
in 152 Dy. While at Q, x = 0.0 MeV the quadrupole moments calculated with and without NM differ 
by ~ 10 %, this difference reaches 2.7% at Q x = 1.00 MeV. The notation of the lines is given 
in the figure. The single-particle orbitals are given from the bottom to the top according to their 
order within the N = 1 shell, see Fig. ||. They are labeled by means of the asymptotic quantum 
numbers [iVn 2 A]r2 (Nilsson quantum numbers) of their dominant component of the wave function 
at £l x = 0.0 MeV and the sign of the signature quantum number r = ±i given as superscript. Since 
the absolute values of (j x )i are small for the [101]3/2 + orbital, the vertical scale of panel (c) is 
different from the other panels. In some panels, the horizontal thin lines with (j x ) = are shown in 
order to guide eye. Top panel shows the proton and neutron angular momenta J2i < 3x{N = 1) >% 
generated by the orbitals of the N = 1 shell. The results of the calculations with and without 
NM are displayed. Note that in the lowest SD configuration of 152 Dy the total J values increase 
smoothly from at n x = 0.0 MeV to 88.26ft (NM) and 70.30?i (WNM) at n x = 1.0 MeV. 
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FIG. 2. The contributions to the expectation values of the spin and orbital angular momenta 
of the particles in the single-particle orbitals forming the N = 1 shell caused by NM. 
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FIG. 3. Neutron single-particle energies (routhians) in the self-consistent rotating potential, 
calculated with and without nuclear magnetism, as a function of the rotational frequency Q x . 
They are given along the deformation path of the lowest SD configuration in 152 Dy. Only all the 
N = 1 and the lowest N = 3 orbitals are shown. The right panel shows the [101]3/2 and [101] 1/2 
orbitals in more details. The notation of the lines is given in figure. 
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FIG. 4. The effective alignments (in h) between the SD configurations in 151 Tb and 152 Dy 
(panel (a)) and in 152 Dy and 153 Dy (panel (b)) are shown by the lines without symbols. The 
configurations in 151 Tb and in 153 Dy are obtained with respect to the lowest SD configuration in 
152 Dy by emptying of the 7r[651]3/2(r = +i) orbital and by the occupation of the 7r[402]5/2(r = +i) 
orbital, see Ref. for details. The expectation values of the single-particle angular momenta 
(jx)i of these orbitals are shown by the lines with symbols. Solid and open symbols are used for 
values calculated in 152 Dy and the neighboring odd nucleus, respectively. Solid and dashed lines 
are used for the results obtained with and without NM. In panel (b), the (j x )i values calculated in 
153 Dy are not shown since they differ from the ones obtained in 152 Dy by less than 0.026?i. 
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